#create matrix
Id <- matrix(data= c( 1.00, 	 0.56, 	 0.55, 	 0.33, 	 0.38, 	 0.48, 	 0.53, 	 0.40, 	 0.42, 
         0.56, 	 1.00, 	 0.48, 	 0.35, 	 0.39, 	 0.52, 	 0.57, 	 0.40, 	 0.41, 
         0.55, 	 0.48, 	 1.00, 	 0.30, 	 0.37, 	 0.44, 	 0.47, 	 0.42, 	 0.39, 
         0.33, 	 0.35, 	 0.30, 	 1.00, 	 0.40, 	 0.32, 	 0.35, 	 0.38, 	 0.36, 
         0.38, 	 0.39, 	 0.37, 	 0.40, 	 1.00, 	 0.37, 	 0.40, 	 0.40, 	 0.39, 
         0.48, 	 0.52, 	 0.44, 	 0.32, 	 0.37, 	 1.00, 	 0.62, 	 0.37, 	 0.38, 
         0.53, 	 0.57, 	 0.47, 	 0.35, 	 0.40, 	 0.62, 	 1.00, 	 0.41, 	 0.41, 
         0.40, 	 0.40, 	 0.42, 	 0.38, 	 0.40, 	 0.37, 	 0.41, 	 1.00, 	 0.56, 
         0.42, 	 0.41, 	 0.38, 	 0.36, 	 0.39, 	 0.38, 	 0.41, 	 0.56, 	 1.00 ), 
  nrow=9,             
  ncol=9,               
  byrow = TRUE)

dimnames(Id) = list( c("Hsap1",	"Ggal1",	"Drer1", "Hsap1L1",	"Drer1L1",  	"Hsap1L2",	"Ggal1L2",  	"Hsap1L3",   	"Ggal1L3"),
                     c("Hsap1",	"Ggal1",	"Drer1", "Hsap1L1",	"Drer1L1",  	"Hsap1L2",	"Ggal1L2",  	"Hsap1L3",   	"Ggal1L3")) # add row and column names

#create heatmap with pheatmap
library(pheatmap)

pheatmap(Id, display_numbers = T,color = colorRampPalette(c('lightcyan','purple'))(100),fontsize_number = 14,fontsize_row=16, fontsize_col=16,border_color=NA,cutree_rows = 2, filename = "Id_DNase1_HM.pdf")

